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fV] ■ ABSTRACT 

• , We present Herschel observations of the water molecule in the massive star-forming region W3 IRS5. The o-H2'0 1 lo- loi , p-Hj'^O 1 1 1 -Ooo, P-H2O 

1-*^ ■ 2o2-lii, P-H2O In-Ooo, 0-H2O 22i-2|2, and 0-H2O 2i2-loi lines, covering a frequency range from 552 up to 1669 GHz, have been detected at high 

Mh' spectral resolution with HIFI. The water lines in W3 IRS5 show well-defined high-velocity wings that indicate a clear contribution by outflows. 

JL [ Moreover, the systematically blue-shifted absorption in the H2O lines suggests expansion, presumably driven by the outflow. No infall signatures 

»^ , are detected. The P-H2O lu-Ooo and 0-H2O 2i2-loi lines show absorption from the cold material {T ~ 10 K) in which the high-mass protostellar 

envelope is embedded. One-dimensional radiative transfer models are used to estimate water abundances and to further study the kinematics of the 

region. We show that the emission in the rare isotopologues comes directly from the inner parts of the envelope {T > 100 K) where water ices in 

the dust mantles evaporate and the gas-phase abundance increases. The resulting jump in the water abundance (with a constant inner abundance of 

, lO""*) is needed to reproduce the o-H,'0 lio-loi and p-H'*0 In -Ooo spectra in our models. We estimate water abundances of 10"** to 10"' in the 

^Nj ■ outer parts of the envelope {T < 100 K). The possibility of two protostellar objects contributing to the emission is discussed. 

> ■ 
1 ^ Key words. ISM: molecules - ISM: abundances - Stars: formation - Stars: protostars ~ Stars: early-type - Line: profiles 

00 : 

C^ ; 1. Introduction overall abundances of less than 10"^ jBoonman et al.l |2003|) . 

r — , _, , , • , ■ r , • ,- ■ The much smaller beam of the Herschel Space Observatory 
••_ he water mo ecu e IS a kev snecies tor stiinvinp star formation . ~„,„, . 



C/2 



, The water molecule IS a key species for studying star formation, (p^bratt et al. 2010), together with the high spectral resolu- 

g . In the cool regions of molecul^ clouds {T < 100 K), water is tion of the Heterodvn e Instrument for the Far-Infrared (HM, 

' present as ice in the mandes of dust grains. In the immediate |de Graauw et al.ll201(ih . gives a unique opportunity to separate 

surroundings of high-mass protostars, the dust is heated to tem- ^^^ ^^^^.^^ ^^^ absorption in the inner (r > 100 K) and outer 

_ _ ^ peratures well above 100 K, evaporating the water ices and in- (7^ < jqo K) protostellar envelope in order to build a model of 

^ , creasmg its abtindance in the gas phase by several orders of mag- ^^^ ^^^^^^^^^ ^^^ kinematics of these innermost regions which 

^ . nitude and making water one of the most abundant molecules. ^^^ ^^.^^^^^j ^^ ^^^j^^^ ^j^^ ^^^^^^ ^^^ 

rrt ' Understanding how the accretion of matter overcomes radiative t . • , ,-,■.• 

. . ^ ^ . • . u ■ 1 ui In this letter, we report water observations of the region 

pressure in massive protos tars is a maior astrophvsical problem r-- — ..p,,. *^ mrTT^;^ nn tdc< • u- u, 

!lZinnecker& Yorkg llo^. and water mav plav a n active role ^^ IRS5 (|Wynn-WdIiams et alj [1971. W3 IRS5 is a bright 

in the energy balance (e.g. IPotv & Netifelil997l) . Due to the !"frt''^^ '"'^'"'^^ ^°'^^'^'^ '?. '^^ active star-for ming region W3 

, , . , J r — f £11 » J in the Perseus arm at a distance of 2.0 kpc Machisuka et al. 



abundance jump, water makes it possible to specifically study ^„„7t, , , . , „.„ , ,, ' ^^-iT — 1 , . . 

,. . • V u- u,u ■ . t \ 2006b. Its high FIR luminosity (1 - 2 x 10^ L^. Laddetal 



the inner regions from which the massive protostar accretes. tt^t^' , , , ,- .... rrrr, — "' . l U^^^J 
, , , .. . u »■ r » » 11 |1993) and the detection of radio emission (^Wilson et al. 2003 
Low angular-resolution water observations of protostellar ' ^ " ' 



1 , ■ ^ ni • ■ -.u J V /IT \ van d er Tak et al. 2005) indicates that W3 IRS5 hosts high-mass 

envelopes (region < 0.1 pc m size with a density n(H9) ' . , \- , • ,.- . , • tt, • 

^ ia6 -3 ■ .11 1 J J ,. J . . stars in an early stage of evolution. High-resolution near-IR im- 

> 10 cm and centrally peaked density and temperature , , ,-t., • t-,,^ tt^o^ • . 

ci N vu ^u J e J c /^u » /TC/-VN J *u ages show a cluster of IR-sources in W3 IRS5 with two sources 

pronles) with the Infrared Space Observatory (ISO) and the . r ■ r- , . ■ , /r.- m , rrr ; n 

e u ir ^ \\r a » c . ii-» /c\jf a c\ u i a identified as the mam elements (Figure 1 and Megeath et all 

SubmdhmeterWave Astronomy Satellite (S WAS) have revealed u^^^ i ^, , , , • , -,7^ 77777 

I2005D . These correspond to the bright millimeter sources MMl 

* Herschel is an ESA space observatory with science instruments pro- and MM2 identified by Rodon et al. ( 2008 ). They are separated 

vided by European-led Principal Investigator consortia and with impor- by ~ I" and drive tw o of the several outflo ws identified by 

tant participation from NASA. iBoonman et all (l2003h and lRodon et all (l2008h . 
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Fig. 1. NICMO S color-composi t e imag e of W3 IRS5 and its sur- 
roundings from iMegeath et all (l2005h . The white circle repre- 



sents the HIFI beam size for the P-H2O In-Ooo and p-Hj^O ln- 
Ooo lines (20")- The gray-scale inset is a close-up of the central 
area at 2.2 microns with white contour s corresponding to 3 0, 60 
and 90% of the peak flux at 1.4 mm (iRodon et al.ll2008h . The 
black bar (top right of inset) represents a distance of 5.000 AU 
(2.5"). 



The first observations of the P-H2O 2o2-lii 
and 0-H2O 2i2-loi water lines and the rare o-Hj^O lio 
p-HfO 



O-H2O 221-2i2 

loi and 
lii-Ooo isotopologues towards a massive protostar are 
presented here. 



2. Observations and data reduction 

W3 IRS5 was observed with HIFI as part of the Science 
Demonstration Phase, between March 3-22, 2010. The position 
observed is at R. A. 02:25:40.6 and Dec +62:05:51 (J2000), clos e 
to the peak radio and mid-IR emission dvan der Tak et al.ll2005h . 
The observations are part of the Priority Science Program (PSP) 
of the Guaranteed-Time Key Program Water In Star-forming re- 
gions with Herschel (WISH; van Dishoeck et al, in prep). 

Data were taken simultaneously in H and V polariza- 
tions using both the acousto-optical Wide-Band Spectrometer 
(WBS) with 1.1 MHz resolution and the correlator-based High- 
Resolution Spectrometer (HRS) with 480 kHz resolution (reso- 
lutions of 0.30 and 0.13 km s ', respectively, for the 1100 GHz 
lines). We used the double beam switch observing mode with a 
throw of 3'. HIFI receivers are double sideband with a sideband 
ratio close to unity. 

The frequencies, energy of the upper levels, system tem- 
peratures, integration times, rms noise level at a given spectral 
resolution, the beam size and efficiency Bgff, and the observed 
continuum level for each of the lines are provided in Table [T] 
Typical uncertainty i n the frequencies is on the order of 100 kHz 
(iPearson et al.lll99ll) . Calibration of the raw data onto T^ scale 
was performed by the in-orbit system (lRoelfsemall2010h : con- 
version to rjnb was done with a beam efficiency estimated by 
raster maps of Mars (Table 1, R. Moreno, priv. comm.). The for- 
ward efficiency is 0.96. Currently, the flux scale accuracy is es- 
timated to be between 5 and 10%. Data was calibr ated in the 
Herschel Interactive Processing Environment (HIPE, IOttll2010h 



version 3.0. Further analysis was done within the CLASSQpack- 
age. After inspection, data from the two polarizations were aver- 
aged. 



3. Results and analysis 

Our interpretation of the spectra uses the following terms: inner 
and outer envelope (as defined in §[TJ, cold molecular cloud (the 
region with T ~ 10 K in which the protostellar envelope is em- 
bedded), and foreground cloud (unrelated molecular material in 
the line of sight at a different Vlsr)- 

The continuum subtracted water spectra are shown in Fig.|2] 
For the o-H'^O lio-loi, P-H2O 2o2-lii, p-H'^O ln-Ooo, and p- 
H2O 1 ii-Ooo lines, we show the HRS spectra. For the 0-H2O 221- 
2i2 and 0-H2O 2i2-loi lines, we show the WBS spectra because 
the velocity range covered by the HRS for those lines is not suf- 
ficient. The 0-H2O 221 -2i2 line is contaminated by a CH line 

(1661.113056 GHz) seen in absorption at 58 km s '. The p- 

H2O lii-Ooo and 0-H2O 2i2-l()i lines show a narrow absorption 
feature at -20.4 km s"'. Since W3 IRS5 is located very close to 
the galactic plane (b ~ 1.2°), this water absorption is probably 
produced by a cold (T ~ 10 K) foreground cloud on the line 
of sight. Other cases of absorption due to cloud s along the line 
of sigh t to ward similar sources ar e reported by Marseille et al.l 
dloTol) and lWvrowski"eran JIoToh . 



3.1. Emission and absorption features 

Most of the water emission spectra in W3 IRS5 can be described 
as the sum of two velocity components: one broad (FWHM 
= 33-40 km s^') and one medium (FWHM = 5-10 km s^'). 
The broad component is visible in all but the rare species and 
is centered approximately on the Vlsr of the source (see the 
broad component for the 988 GHz line in Fig.l2Ti. lBoonman et al.l 
(120031) observed the outflow in W3 IRS5 in CO 7 = 7 - 6 with 
the JCMT and found a similar line width to what is in the water 
broad component. We assume that the broad component is due 
to the outflow. 

The medium component exhibits emission in the P-H2O 
2o2-lii and rare isotopologue lines with two peaks at -37 and 
-41 km s"' in the P-H2O 2o2-l 11 and p-Hj^O 1 n-Ooo spectra. A 
blend of emission and absorption is seen in the P-H2O 1 1 1 -Ooo, o- 
H2O221-212, ando-H20 2i2-loi lines. The absorption in the H2O 
lines is blue-shifted (P-Cygni profi le), suggest i ng an expansion 
of the envelope, as also detected bv lBenz et al.l (1201 Ol) in the hy- 
dride lines of W3 IRS5. The expansion is probabl y powered by 
the multiple outflows known to exist in W3 IRS5 (iRodon et al.l 
2008). The P-H2O ln-Ooo and 0-H2O 2i2-loi lines have a dou- 
ble plateau in absorption at -38 and -42 km s"'. The absorption 
at -38 km s~' is highly saturated. A simulation using RADEX 
(Ivan der Tak et al .120071) for a cold cloud with a kinetic tempera- 
ture of 10 K, A^HzO = 3 X 10'^ cm-2, a fine FWHM of 2.8 km s"', 
and «H2 = 1 X 10^ cm"-^ gives opacities of 4.0, 2.2, 2 x 10"^ 
and 3 x 10"^, respectively, for the P-H2O ln-Ooo, 0-H2O 2i2- 
loi, 0-H2O 22i-2i2, andp-H20 2o2-lii species. This agrees with 
the detection of saturated absorption only in P-H2O 1 1 1 -Ooo and 
0-H2O 2i2-loi, as well as the absorption by a cold foreground 
cloud at -20.4 km s ' also detected only in those lines. The opti- 
cal depth of the other water lines in cold regions is so low that no 
absorption is visible. The absorption at -38 km s"' presumably 
comes from the cold molecular cloud associated with the source. 



http://www.iram.fr/IRAMFR/GILDAS 
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Table 1. Lines of H2O observed with Herschel/HIFI in W3 IRS5. 



Water species 



Frequency 
(GHz) 



(K) 



i sys 
(K) 



fin, 6v rms Beam Vrnb cont. 

(s) (kHz) (K) (") (K) 



o-H^'Olio-loi 


552.020960 


61.0 


73 


158 


250 


0.06 


40.0 


0.72 


0.35 


P-H2O202-IU 


987.926764 


100.8 


347 


386 


500 


0.07 


23.0 


0.71 


1.4 


p-Hf l„-Ooo 


1101.698256 


53.4 


351 


601 


500 


0.09 


20.0 


0.71 


2.15 


P-H2OI11-O00 


1113.342964 


53.4 


351 


601 


500 


0.09 


20.0 


0.71 


2.15 


0-H2O221-21, 


1661.007637 


194.1 


1609 


1164 


250 


0.25 


13.5 


0.69 


5.5 


0-H2O212-I01 


1669.904775 


114.4 


1609 


1164 


250 


0.25 


13.5 


0.69 


5.5 




ViisH [km/s] 



,p-H^«0 1„-Ooo 



Fig. 2. HIFI spectra of water lines observed in W3 IRS5. Columns 1 and 2 show the spectra for the o-Hj^O 1 lo-lo 
andp-H20 2o2-lii lines, columns 3 and 4 show the spectra for the P-H2O ln-Ooo, 0-H2O221-212 ando-H20 2i2-loi lines. The single 
protostellar envelope model fits are shown as red lines over the spectra in columns 1 and 3. The two protostellar envelopes model is 
shown as red lines over the spectra in columns 2 and 4. The dashed red line over the P-H2O ln-Ooo ando-H20 2i2-loi spectra show 
the corresponding model without the cold molecular cloud (see § 13. 11 1. The dashed black line over the o-H^^O 1 lo-loi and p-Hj^O 
ln-Ooo spectra shows the model without the abundance jump. Vertical dotted lines indicate the model Vlsr at -38 km s"' for the 
single protostellar envelope model (columns 1 and 3) and -37 and -40.5 km s ' for the two protostellar envelopes model (columns 
2 and 4). Vertical dashed lines in the 0-H2O 221 -2i2 spectra show the position of the CH line. As an example, the outflow component 
for the P-H2O 2o2-lii line is shown in blue. The temperature scales are shown on the left side of the figure for columns 1 an 2 and 
on the right side for columns 3 and 4. 



The two-component profiles seen in emission in several of 
the water lines seem to be systematic features of protostellar ob- 
je cts, since it has been ob served for low-mass Class objects 
bv lKristenseii et all (l2010h and for intermedia te-mass protostars 
bv I Johnstond (|20 1 Ol) . iKristensen et al.l (1201 Ol) propose that both 
the broad and medium components are associated with outflow 
shocks. In contrast to the low-mass protostars, the emission of 
the rare isotopologues only comes from the medium component 
and not from the broad outflow. 



3.2. Two protostellar objects seen in the optically thin lines? 

Double peaks in emission at Vlsr ~ -36 and -40 km s"' in 
the mostly optically thin P-H2O 2o2- 1 1 1 and p-H'^O 1 1 1 -Ooo line 
profiles suggest the detection of more than one pr otostellar en- 
velope . Since the interferometric observations from lRodon et alj 
(l2008h show that there are two millimeter sources within the 
HIFI beam (Fig. [T]) that dominate the dust continuum emission 
in this region, it is possible that the two peaks correspond to 
emission from these two high-mass protostellar objects. 

Ground-based obser vations of C^^O , CS, SO, and 
CH3OH molecules (Ivan der Tak et al.1 I2003L 120061: 



iHelmich & van Dishoeckl Il997h show lines that cover the 
entire velocity range of emission, but they do not help dis- 
tinguish between one or t wo protostel la r env elopes. The 
continuum observations by iRodon et al.l (l2008h could not 
measure the velocities of MMl and MM2. 

Even though the hypothesis of the detection of two proto- 
stellar envelopes in the spectra needs to be investigated in more 
detail (to be presented in a forthcoming paper), we show in § 13.31 
that a two protostellar envelope model tends to fit the observa- 
tions better than does a single protostellar envelope model. 



3.3. A model of W3 IRS5 



We used the Monte Carlo code MC 3D (IWoljl2003l) and the ra- 
diativ e transfer program RATRAN dHoeerheiide & van der Tala 
l2000h to model the water emission in W3 IRS5 following the 
method described in [Marseille et al.l (|2008|) with a power-law 
density exponent of -1.2. The models assume spherical symme- 
try. A more detailed investigation of the region's morphology 
using 2D models is in progress and will be published when the 
entire set of water lines (including maps) is observed. 
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Table 2. Parameters of single and two protostellar envelope (PE) 
models. 



Parameter 


Single PE 


PEl 


PE2 


Outer radius (AU) 


12000 


12000 


12000 


Luminosity (L©) 


10^ 


5x10^* 


5x 10* 


Mass (Mq) 


250 


125 


125 


%20 


2 X 10-** 


1.8 X 10-** 


8 X 10-'" 


Post-jump Xjj^o 


1 X lO-t 


1 X 10-* 


1 X 10-* 


VturCkms-') 


2.0 


0.4 


0.5 


Vexp (km s"') 


2.0 


1.1 


1.2 


yLSR(kms-') 


-38 


-37.0 


^0.5 



Our models in W3 IRS5 have three components: an out- 
flow, one (or two) protostellar envelopes, and a cold cloud. The 
outflow contribution to the spectra was modeled as a Gaussian 
emission line with FWHM between 33 and 40 km s"' . The cold 
cloud was introduced taking the cloud parameters described in 
§ 13.11 for the saturated absorption. Three input parameters are 
used to fit the line profiles in the envelope model: water abun- 
dance (Xn^o)^ turbulent velocity (Vtur), and expansion velocity 
(Vexp)- The width of the line is adjusted by changing Vtur- The 
line asymmetry is reproduced by the expansion velocity Vexp- 
The line intensity is fit by a combination of the abundance, Vtur, 
and Vexp parameters. The same Vtur and Vexp have been assumed 
for all lines. We use abundance ratios of 500 for H2O/HPO, 



4.13 for H^^O/H^^O (as measured in the W3 region by iPenziasI 
1 198 It IWilson&koodlll994t IWouterloot et"aDl2008h . and 3 for 
ortho/para-H20. The models assume a jump in the abundance 
in the inner envelope at 100 K with a constant inner abundance 
of 10 "*. Having the jump occur at the radius predicted by the 
physical model leads to overestimating the water line intensity 
in the inner envelope. We thus decreased the size of the emitting 
region (to ~1500 AU) by putting the abundance jump at a tem- 
perature of 150 K. An alternative model leading to similar results 
keeps the jump at 100 K and decreases the water abundance in 
the inner envelope by a factor of 10. We assume that the line 
intensities from the various components can simply be added as 
in the optically thin limit, and future models will treat this effect 
in more detail. Table |2] gives the parameters used in the models 
and the model spectra are shown in Fig.|2]for both a single and a 
double protostellar envelope model. The outer radius in the table 
corresponds to the deconvolved radius taken from the 1.3 mm 
continuum emission (Oldham etal. 199 4). 

The observed lines were fit using the following strategy. 
First, we modeled the rare isotopologue lines and then the P-H2O 
2o2-lii line (including the outflow emission) using the same 
^H^O. ^tur, and Vexp values. Once we were able to reproduce 
the main features of those profiles by minimizing residuals on a 
grid of values, we applied the same parameters to the rest of the 
lines, including an outflow component when justified. For the 
single protostellar envelope model, it was not possible to cor- 
rectly fit the rare isotopologues and P-H2O 2o2-lii lines using 
the same Xh,Oi ^tur, and Vexp values. For the double protostellar 
envelope model, all of the spectra are better fit with the excep- 
tion of 0-H2O 221 -2 12. The deep absorption in the P-H2O 1 n-Ooo 
and 0-H2O 2i2-loi lines is reproduced well only if we add a cold 
cloud to the model. The uncertainties for the water abundances 
are about a factor of 3 and for Vtur and Vexp about ~ 1 km s"'. 

As shown in Fig. 2, the emission in the rare isotopologue 
lines is only reproduced in our models by including a jump in 
the water abundance in the inner envelope. The estimated op- 
tical depths for the lines of the rare species are low (r ~ 0.1, 
estimated using RATRAN). This suggests that the emission in 



the rare species comes from the inner envelope where the wa- 
ter abundance is greatly enhanced. Emission from this region 
could be the expected contribution from the passively (radia- 
tively) heated inner envelope. 

4. Summary 

Spectra of six water lines observed with the Herschel-HIFI in- 
strument in the high-mass star-forming region W3 IRS5 are pre- 
sented and discussed. AID radiative transfer model was applied 
to estimate water abundances and study the kinematics in the 
protostellar envelope. The main results follow. 

1 . We detect strong water emission and absorption. 

2. A strong outflow component is detected in the water lines. 

3. Water emission shows absorption from cool molecular gas 
hosting the protostellar envelope. Blueshifted absorption 
suggests an expansion of the outer envelope, and no infall 
signature is detected. 

4. Radiative transfer models indicate water abundances ranging 
from 10"** to 10"'" in the outer envelope. 

5. Based on our model, a jump in water abundance in the inner 
envelope is needed to reproduce the H^^O and Hl^^O lines. 

6. The optically thin line profiles are better fit using a model 
with two protostellar envelopes, in agreement with previous 
interferometric continuum observations of W3 IRS5. 
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